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We report on the theory of logarithmic temperature profiles in very strongly devel- 
oped thermal convection in the geometry of a Rayleigh-Benard cell with aspect ratio 
r = 1 and discuss the degree of agreement with the recently measured profiles in the 
ultimate state of very large Rayleigh number flow. The parameters of the log-profile 
are calculated and compared with the measure ones. Their physical interpretation 
as well as their dependence on the radial position are discussed. 
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I. INTRODUCTION 



The ultimate state of thermal convection in a Rayleigh-Benard cell is generally believed 
to be fully turbulent, both in the bulk and in the boundary layers and with respect to both 
the velocity as well as to the temperature fields. As shown in ref. this implies very 
characteristic scaling behavior of the heat transport flux, Nu oc Ra ' 38 as well as of the 
strength of the large scale flow, Re oc Ra ' 50 . Here the Rayleigh number Ra, the Nusselt 



number Nu, and the Reynolds number Re are defined as usual, see e.g. 



of recent experiments, cf. 



3j. In a series 



4-8], these scaling laws of Nu and Re have been measured and 
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reported. 

In the previous work 1[ we have calculated the global scaling exponents in the ultimate 
state. This ultimate state theory, though derived by employing the characteristic profiles of 
fully developed turbulent flows, concentrated on the implications for and the interpretation 
of the various global scaling exponents. We did not explicitly describe and report the local 
profiles themselves. Meanwhile the local thermal profiles in the very large i?a-regime have 



been measured, cf. 



9]. We thus present here the corresponding local profiles, in the spirit 



of and in extension of our earlier theory from ref. 1] 



II. THERMAL LOG-PROFILES OF FLOW ALONG PLATES 

We start from the equations of motion for the thermal field T(x, t). The involved velocity 
field is understood to solve the corresponding Navier-Stokes equation for an incompressible 
flow, V • u — 0, i.e., we assume that the temperature behaves like a passive scalar. The 
molecular properties kinematic viscosity v and thermal diffusivity k are taken as temperature 
independent fluid parameters. It is 

d t T = -u ■ V T + kAT, (1) 

together with the respective boundary conditions (b.c.) T(z = 0) = and T(z = L) = —A 
and no heat flux at the sidewalls. In order to understand its turbulent solution we model 
the RB flow in a rather simplified way, which in the first steps allows to make some use of 



textbook wisdom, cf. jlfjf| §54, or [11] §7.1.4. Consider a flow along an infinite plate in re- 
direction and ask for the profile in z-direction perpendicular to the plate. First decompose 
the fields into their long time means U, O and their fluctuating parts u', 8'. The time 
averages depend on z only and the velocity has an x-component only, choosen to be in 
positive x-direction. Per definition we choose 0(z = 0) = 0. If there is heat flow upwards, 
J = (U3O) A,t — 1^83(8) A,t positive, 9 decreases with increasing distance z from the plate and 
vice versa. This describes the ! situation at the bottom plate of an RB-cell. - At the top 



plate the temperature flux J is still positive, directed towards the plate, but 9 increases 
with growing distance from the top plate or decreasing z, starting from z = L downwards. 



o = -u ■ v e - u' ■ w +kAq. (2) 

The first term vanishes, since only the x-derivative can contribute but 6 depends on z only. 
The last term for the same reason reads nd^Q(z). The second term has to be modeled; we 



use the mixing length type Ansatz for an eddy thermal diffusivity. First, u'fi' as a long time 
mean depends on z only, is approximated by employing the temperature gradient hypothesis, 



u'fi 1 = — Kturb{z) d z Q(z) for i = 3, otherwise . (3) 

The eddy diffusivity K turh (z) is supposed to be a property of the flow, not of the fluid. The 
physical flow properties on which K tur b(z) can depend, are the distance z from the plate and 
the characteristic velocity scale u* of the turbulent fluctuations, defined by the wall stress 
vU x \ z (z — 0) = ul, the shear rate or drag on the wall. From dimensional reasons we write 

Kturb=K>8 Z . (4) 

Here the factor R g is the dimensionless thermal von Karman constant, whose empirical value 
depends on the flow type and for Rayleigh-Benard flow is not known. 

The equation for the temperature profile then is d z [0+K turb (z)d z <d(z)+K,d z <d(z)]. Integrate 
from z = to some arbitrary z and use n tU rb{z = 0) = to find (ftturb + n)d z Q(z) — 
nd z Q(z) |o = 0. The last term is just J or Nu ■ /cAL -1 , where A = T b — T t is the temperature 
difference between the bottom and the top plates, which causes the thermal flow. J is 
^-independent. This leads to the thermal profile equation 

d-Mz) = — . (5) 

K + Kq Z -U* 

If J is positive, i. e., the heat flows upwards, the thermal gradient is negative, Q(z) decreases, 
as is characteristic for the bottom plate. With the choice Q(z = 0) = at the bottom plate 
we have Q(z) = T(z) — T b relating 6 with the physical temperature T(z). 



From this equation we draw conclusions for the bottom part profile. In the immediate 
vicinity of the plate it is k ^> Kturb, giving rise to the "linear thermal sublayer" of the profile, 
Q(z) = — J ■ (z/k), < z < z* jK . This linear thermal sublayer extends until z = z* )K for 
which k=kq z* 7K , i.e., if the molecular and the eddy thermal diffusivity are of the same 
size. The relation to the kinetic sublayer width is 

z* K = — — -)■ z* K = Pr^ 1 ^— = Pr^ 1 — , with z* = — . (6) 

KgU* KgU* Kg 

Beyond the linear thermal sublayer the profile is increasingly dominated by the eddy 
diffusivity, K turb ^> k, implying d z Q(z) = —J/(R e u*z) and thus 

Q(z) = -— fin- + f) (7) 

is the thermal profile in the turbulent BL. The integration constant / may depend, of course, 
on k and z/, but only in the form u/k, since / has unit 1. Thus / = f(Pr), which is reported 



m |10| to empirically be about 1.5 for air, i. e., for Pr pa 0.7. 

To make use of this formula for insight into the thermal log-profile we need the strength 
of the velocity fluctuations function of Re(Ra). This has been derived in But 

before determining u* we make now contact to the experimentally measured profile, cf. 9]. 
The experimental log-profile is parametrized in the form 

T ff ~J m = A In 4 + B . (8) 
T b -T t L 1 ; 

Here T m = (TJ, + T t )/2 = T b — A/2 is the mean temperature in the Rayleigh-Benard cell. 
Comparing eqs.flTJ) and (jHJ) leads to 

We now can identify the dimensionless empirical parameters A and B. 

kNu , , 

A = -— — 10 

LKgU* 



and 



1 kNu 
B — — 

2 Ln e u* 



In— + f(Pr) 



A 



In— + f{Pr) 



1 

+ 2' 



(11) 



One easily verifies that the dimensions of A and B are indeed 1 . Besides trivial parameters 
and an yet undetermined empirical parameter, the thermal von Karman constant fig, two 
physical quantities determine the measured constants A and B. This is first the strength 
of the heat flux J = Nu ■ kAL -1 , describing the amplitude of the log-profile, and second 
the strength of the turbulent velocity fluctuations it*. While Nu and its scaling behavior 



in the ultimate range is experimenta 
amplitude has been calculated in 



ly (and theoretically) rather well known, the velocity 
l| . We make use of those results now. 



III. THE LOG-PROFILE PARAMETERS A AND B 
A. Parameter A 

We start with the discussion of the amplitude A of the log-profile. From eq.f llOp the 
parameter A apparently can be written as 

kNu Re^Pr^Nu 

LRgU^ KgU^/U 

with Re = UL/v. Here we have introduced the wind amplitude U and the corresponding 
Reynolds number Re. For not yet too large Ra this amplitude U usually is visualized as 
a large scale circulation (LSC) in the RB sample. For very large Ra on the other hand, 
which are considered here, such LSC will probably not survive under the strong turbulent 
fluctuations. But its remnants locally in space and time still must have physical importance, 
since there apparently is enough shear in the plates' boundary layers to induce transition 
to turbulence, as it can be observed experimentally in the scaling exponents of the Nusselt 
number Nu and the Reynolds number Re versus Ra as well as the measured characteristic 
changes of the scaling exponents indicating this transition, see |4do| and also [9]. 
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To quantify this we use results from 



121 ] . In the ultimate range Ra > 10 15 the effective 

39 x 



Reynolds number was found to be Re = 0.0439 xRa ' 50 , resulting in Re(Ra = 10 15 ) = 
10 6 . Extrapolating Re from smaller Ra, known as the classical range of RB flow, in 



12| it 



is found Re = 0.407i?a a423 , the scaling exponent being well consistent with the GL theory 



131 ] . According to this classical range formula the wind amplitude would be measurably 
smaller at Ra = 10 15 , namely Re = 0.901 x 10 6 . 

In the cited RB experiment it is Nu = 5631 and Pr = 0.859 at Ra = 1.075 x 10 15 . 
From this we can calculate the coherence length in the turbulent bulk, which is defined as 
i coh « l0r) Kol = lO^/e) 1 / 4 = \0(Pr- 2 RaNu)-- l l A L. One gets £ coh /L = 1.87 x 10" 4 or in 
physical units (using L = 2.24 m for the height of the Goettingen Uboot device) £ co h = 0.419 
mm. This sets the lower bound of the turbulence eddy sizes in the bulk. There will be larger 
eddies too, whose extension is between the external scale of order L (= D in the case T = 1) 
and l co h- These even if fluctuating temporally and in position can provide the necessary 
shear in the boundary layers. 

To quantify this we compare with the thickness of the linear viscous sublayer above 
the plate, which will be introduced and calculated later. It will be estimated as z*/L = 
1.98 x 10~ 5 for Ra = 10 15 . Then z*/£ co h = 0.106, meaning that even the smallest eddies 
of the bulk cascade are still 10 times larger than the kinetic viscous sublayer extension. 
That holds all the more for the energy carrying larger eddies up to the macroscocic scale L. 
Thus the bulk flow, even if strongly fluctuating, can provide sufficient shear to imply the 
turbulence transition at the plates (and most probably also at the side walls). 

Having clarified the meaning and physical importance of the wind amplitude U, we can 
determine the log-profile and its parameters A and B. In [I] we have derived the following 
two expressions for the Nusselt number Nu and the velocity fluctuation amplitude relative 
to the given asymptotic flow velocity, u*/U, in the logarithmic ultimate range, cf. eqs. (20) 



and (2) of that paper, namely 



^RePr , 

Nu = - 2 g- = . 13 

In (tfef \) +f(Pr) 

The fluctuation amplitude solves the equation 

^ = (u) 

Here h is an empirical parameter of the velocity profile. It characterizes the position of the 
buffer range, in which the transition occurs from the linear viscous sublayer to the log-layer 
range in the velocity profile. Usually that profile - known as the law of the wall - is written 
in the form 

um i ln w +5u _ (15) 



Here R is the kinetic von Karman constant, taken in the following as R = 0.4. B u describes 
the velocity U(z) at z — z*, i. e., it characterizes the buffer range of the velocity profile 
by indicating, how far the log-law velocity U(z) is shifted in amplitude, if extended to the 
linear viscous sublayer. There are various values of the empirical constant B u given in the 



literature, such as e. g. 5.1, cf. 10], JjJ. To the best of our knowledge there is not yet any 
information available for its value in the case of Rayleigh-Benard flow in a closed container. 
Note that this RB-flow neither is stream wise infinite nor does it span wise approach a given 
constant amplitude U, but rather becomes small again near the mid plane and then even 
changes sign. Thus B u remains a parameter still to be measured in RB cells. We include it 
into the profile equation (I14p in the form 

b = e~ RB - . (16) 



According to eq. ffM]) the numerical value of b determines the amplitude of the velocity 
fluctuation scale (or shear stress oc v%) together with the Reynolds number Re in the 
combination Re/b. From (I14j) we see that the smaller b is the smaller will be. Similarly 
will decrease with increasing Reynolds number Re and thus with growing Rayleigh number 
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Ra. For Re(Ra) we meanwhile have experimental information, see [12|: In the ultimate 
state it is Re e ff — 0.0439-Ra 0,50 . We thus have typical .Re-numbers of order 10 6 or more in 
the ultimate range of thermal convection with Ra ~ 10 15 . 

The velocity profile parameter b also enters the Pr-dependence of the temperature profile 
via the definition f(Pr) = f(Pr) + ln6/2. Using all these previous results we can express 
the amplitude parameter A in the form 

1 1 



.4 



(17) 



2 [in (fRel) + f(Pr)\ 2 [k£ + f(Pr) 

Two remarks may be useful. In this expression for A the thermal von Karman constant Re 
does not appear explicitly; instead the kinetic one, R, shows up. This comes from eliminating 
Nu/Rq in the defining equation ( |T2|) for A by using eq.( lT5|) . Second, the originally derived 
explicit iV-u-dependence can be completely substituted and expressed in terms of the log- 
corrections originating from the log-profile of the velocity. The only reminder to the thermal 
profile is the thermal buffer range shift f(Pr), the thermal analog of the corresponding 
kinetic shift B u . 

Both A and B depend on Re as well as on «*/£/, and thus on Ra. This implies the 
i?a-dependence of these measured fit-amplitudes A and B. Numerical values for Re and 
u*/U have been given in Table I of reference lj] for Ra = 10 14 and 6=1. For any other Ra 
(and b) they can be computed by solving the implicit equation (fl4l) . either numerically or 
via its continued fraction representation given in [l|. We do this below for Ra = 10 15 . 

Taking for the parameters / and b the values reported in textbooks for channel or pipe 
flow of gases, i. e., f(Pr = 0.7) w 1.5 (cf. |l0j) and b = 0(1), one finds f(Pr) = 0(1). The 
estimate for the magnitude of the relative fluctuations u*/U in [l| gave values of the order 
of 0.05. Thus the first term in the denominator of eq. ( JT7|) is an order of magnitude larger, 
even more so for larger Re(Ra). Neglecting therefore / we find the following approximate 



formula for the parameter A: 

A ps J-^ = -J-^. (18) 

21n(ff.Re±) 2k 17 v ' 

For i?a = 10 14 it is u*/U = 0.039 and 0.065 for the choices Re\ and Re<i in table I of l|, 

respectively, leading to the numerical values A = —0.0488 and A = —0.0813, respectively. 

For Ra = 10 15 and Re = 1.4 x 10 6 we calculate from eq. (H4jl (and the choice 6=1) the value 

u*/U = 0.036, leading to A = —0.045, still about half an order of magnitude larger as the 

measured value A exp = 0.0082, given in [9J. 

The scaling of the profile parameters A and B with the Rayleigh number Ra can be 

determined as follows. First transform the .Re-dependence of u+fU into a scaling Ansatz, 

u*/U oc Re a ' . Here a' = a' (Re) is a local exponent, decreasing with Re and therefore with 

Ra. In the relevant .Re-range one calculates a' = —0.0855. The negative scaling exponent 

a' means that the ratio u*/U = Re^/Re decreases with increasing Re and corresponding 

increase of Ra. But, of course, the turbulent fluctuation scale i?e* = u*L/v itself increases 

in size for stronger thermal driving Ra. It is 

.Re, oc Re l+a = Re 9145 , in the i?a-range 10 15 . (19) 

Introducing now this .Re-scaling of i?e* into eq.ljTBl results in 

A oc Re a ' oc Ra°- 50a ' = Ra~ 0M3 , in the ^a-range 10 15 . (20) 

The amplitude A thus decreases with Ra, though very slowly. Over the next decade it 
will be smaller by approximately a factor of 0.906, i.e., its value at Ra is expected to be 
about 9.4 % less. This predicted slow decrease of A with Ra could be consistent with 
experimental observation 9j. We emphasize that a' decreases even further with increasing 
Re, i.e., asymptotially in Ra the amplitude A will approach a constant, i?a-independent 
limit Aoq. 

One may wish to try coming closer to the experimental value for A in the measured range 



of Re ~ 10 6 and Ra = 10 



15 



12j |) by adjusting the still badly known value of b. Instead of 
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b =0(1) one could reduce it, mimicing thus a larger Re, which leads to smaller u*/U . To be 
more precise: Take the measured u*/U = 0.0082 and Re(Ra = 10 15 ) = 1.4 x 10 6 and insert 
this into f fT4|) . This then determines b; we obtain ln(l/6) = RB U w 40, thus B u « 100, and 
b close to 0, namely b = O(10~ 17 ). 

Such unusually large value of B u indicates that the model for the velocity field's [/-profile 
is indeed oversimplified. An important missing feature when assuming an asymptotically 
constant flow amplitude U with increasing distance z from the plate is as follows. The 
LSC in an RB-sample does not approach a nonzero constant value when going off the plate. 
Instead U(z) not only will decrease again, it even passes through at z/L = 0.5 changing 
sign in the upper half of the RB cell. Thus between about 0.25 < z/L < 0.5 the average flow 
has a negative instead of a positive or zero z-slope. The thermal flux J keeps its positive 
sign. Because of this quite different U(z) profile in comparison to pipe or channel flow we 



e in the range above z/L > 0.25. 
9j, Fig.l. - Another consequence 



have to expect a deviation from the temperature log-profile 
Indeed such deviation can be seen in experiment, see 
of such deviation from the p! ure log-profile is that the intimate connection between the 
amplitudes A and B, to be discussed in the next chapter, will loose validity. 

B. Parameter B 

Let us now analyze the parameter B as derived in eq. (TTTj) . in particular the right 
hand part of this equation. Consider first the term ln-^. With the above given defini- 
tion z* = v/u* one obtains L/z* = Lu*/v and thus B = A [in (i?e|f) + f(Pr)] + 1/2 = 



.4 



In (Rejf- 1) + f(Pr) + ln2 +1/2. Now, the ln-term plus the next one are just —jj, which 



after multiplication with A leaves —1/2, canceling the final term 1/2. Thus 

B = A-\n2 or B/A = ln2 = 0.693... . (21) 

Irrespective of any approximate calculation of the coefficient A, the other coefficient B is 
always about 70% of A, in particular is also negative. Using this the T-profile can be written 
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as 



T(z) - T r 
T b -T t 



A (hi- + ln2 



) 



A In 



L/2 ' 



z 



(22) 



This is consistent with the underlying idea that T(z = L/2) = T m , which has been used 
when relating the T-profile, expressed in terms of the thermal current J or Nu, with the 
driving temperature difference A = T& — T t , and which also is an immediate consequence of 
the Ansatz ([8]) for fitting the data. 

Let us come back to the consequences of the deviation of the (^-component of the) velocity 
profile from the classical case with asymptotically constant amplitude U. In an RB cell the 
U(z) -profile does not stay asymptotically constant with z but goes through a maximum, then 
decreases further and even changes sign at z/L ~ 0.5. This different maximum-and-beyond 
behavior of the wind leads to deviations in the thermal profile from a log-law. Therefore 
the intimate relation between A and B, viz. B/A = ln2, typical for the log-profile, will no 
longer be valid. Also this statement is consistent with experimental observation, cf. fsj]: the 
ratio B/A has quite some scatter, and although not too far way it apparently differs from 



To briefly summarize, the thermal profile is characterized by - starting from the plate 
- (i) a tiny linear thermal sublayer of extension an oder of magnitude less than the bulk 
coherence length, (ii) a buffer range in which the linear increase in the sublayer turns over 
into the (iii) log-law, observable over a broader z-range up to about a quarter of the RB cell 
height, then (iv) changing the profile again due to the decrease of the blowing wind including 
its directional change, from positive to negative (or vice versa), which may be denoted as 
the temperature's center profile. This latter one, the center part, still has to be explored in 
more detail. 



In2. 
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IV. POSITION DEPENDENCE OF PROFILE PARAMETERS A(r) AND B[r) 

The basic assumption of the just given derivation of the profile parameters A, B is a 
plane parallel homogeneous flow with veloctity U over an infinitely extended plate. In a 
Rayleigh-Benard cell this -at best- is realized in the center-range of the circular cell, i. e., 
at r = 0. We now make a crude model for the shape of the large scale circulation (LSC). 
We have in mind the case of aspect ratio T = 1, but one can argue similarly for T = 1/2 
(and other). Then, if one moves away from the center range at r = 0, the relevant flow 
velocity leading to the velocity shear in the BL and the perpendicular logarithmic profiles 
of velocity and temperature is the x-component of the LSC only. We therefore have to 
substitute in above formulas always U =>- U x = U x (r). If we consider for simplicity a circular 



LSC, we get U x = Ucoscj) — Uyl — j^-; here U still denotes the LSC amplitude, which 
defines Re = UL/u. 

Considering expression ( fl8l) will result in an r-dependent profile coefficient 



Air) = — — , with A the coefficient at the center. (23) 

Expressed in terms of the (relative) wall distance £ = thus < £ < 1 between the wall 
and the center, respectively, it is 

*® = w m 

The coefficient decreases with increasing distance from the wall oc £ _1,/2 . This result 
is consistent with the experimental finding of an A-decrease towards the center. The data 
presented in [9] have been taken at the position {R — r)/L = 0.0045, corresponding to 
£ = 0.009. This implies A(^ = 0.991) = A(£ = 0.009) = 7A7A. 

We note that these last formulas cease to be valid in the limit £ — > or r — > R, because for 
sufficiently small U x (r) the BL locally is no longer turbulent, finally there is only upward flow 
and one is in the BL of the side wall. But note further that with the full, non- approximate 
expression (fTTl) for the profile coefficient A the additional term in the denominator weakens 
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the r- dependence, all the more the larger r becomes. Then even the limit r — > R exists, 
leading to 

A(r -+R) = --J— . (25) 
2f(Pr) 

One can draw various conclusions from the r-dependence of A(r) on the shape of the LSC 
as follows. Note first that the magnitude of the LSC wind is everywhere the same along its 
closed trajectory, called U or in dimensionless form Re. The direction of the wind varies 
along the LSC orbit. Its local unit vector be called t(s), which is the local tangential unit 
vector along the (closed) LSC curve, s denotes the arc length along the LSC curve. Then 
the velocity relevant for A(r) is given by U x — e x • t U, the x-component of the local LSC 
velocity vector Ut. Given a model for the LSC-curve, e. g. a circle or an ellipse, this can 
be described in terms of a suitable parameter <fi by x(<j>). Expressing the (arbitrarily given) 
parameter in terms of the arc length s, the tangential vector then is the derivative of the 
curve, i.e., t = Measuring on the other hand U x (r) via A(r) all! ows to re-construct 

the LSC curve. 

Some qualitative features of the r-dependence of A(r) are the following: (i) For small r 
the amplitude A(r) always increases as A(r) = A(l + const(-^) 2 ); there is no term linear 
in r because of analyticity reasons. (For a strictly circular flow it is const = 0.5). (ii) For 
r near the side wall, i. e., for small £ = (R — r)/R, the amplitude A{£) ~ £~™ decreases 
with increasing distance from the wall with an exponent n — 1/2 for a circular LSC; for an 
elliptically shaped wind curve it will be steeper, thus n is larger. In the particular case of 
aspect ratio T = 1/2 there may be two circular rolls above each other, which again would 
lead to n — 1/2. In case there is only one single roll this will be elliptically shaped and 
thus n is larger. In case both LSC shapes are present part of the time, the exponent n 
will be somewhere in between, i. e., in any case one would find a steeper decrease of A 
with the ! wall distance as compared to the circular case for Y = 1. This is consistent 
with the measurement of n w 2/3 in the case T = 0.5 (Goettingen Uboot team, private 
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communication), (iii) If T is a little below 1, n will be a bit larger than 1/2; if T is a little 
above 1, n is expected to be somewhat smaller than 1/2; etc. 



V. CONCLUDING REMARKS 

We have calculated the profile parameters A and B, defined in eq.flH]), of the experimen- 
tally measured logarithmic temperature profiles in the ultimate state of Rayleigh-Benard 
convection as realized for very large Ra. In the case of a pure thermal log-law as well as 
reflection symmetry with respect to the middle plane z = L/2 = R oi the T = 1 cylindrical 
sample, which we expect and which should hold for approximately temperature independent 
material properties, the coefficient B equals A up to a factor In 2. Since the real wind profile 
is different from the standard case of approachig a constant when going off the plate, in the 
RB cell instead going down towards the center, even followed by a directional inversion of 
the wind direction, one finds deviations from this value. 

The amplitude A physically measures the strength of the turbulent velocity fluctuation 
scale relative to the LSC velocity U (together with the von Karman constant k); if 
one is sufficiently near to the side wall, it also measures the Prandtl number dependent 
temperature profile shift constant f(Pr). We explain the r-dependence of A(r) by the 
decreasing magnitude of the local flow velocity U x (r) parallel to the bottom plate with 
increasing distance r from the center, which by ( TT8l) or (JTTj) leads to an increase of the 
^-amplitude with increasing r. 

The numerically obtained value for A for the case b = 0(1) does not coincide too well 
with the measured one, if we use the values of the fluctuation scale calculated in [1]. 
When we published that work the meanwhile measured values for the sample's LSC-response 



Re = Re(Ra) (cf. 12j) was not yet known and assumed to be smaller; if our interpretation 
is correct it means, that the experimental u*/U values are smaller in an RB-sample, because 
U oc Re is larger than assumed previously and b is smaller. 
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Thus there is very interesting information in the measured T-profile parameters A(r) and 
B(r). But, as usual in Rayleigh-Benard flow, to complete, check and confirm theoretical 
interpretation and explanation will need information also about the velocity profile and its 
expected considerable deviations from the model of classical channel or pipe flow as is used 
here for a lack of better. 
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